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Abstract 

Background: The evaluation of intestinal trophism, mainly the mucosal layer, is an important issue in various 
conditions associated with injury, atrophy, recovery, and healing of the gut. The aim of the present study was to 
evaluate the kinetics of the proliferation and apoptosis of enterocytes by immunohistochemistry and to assess the 
complexity of intestinal mucosa by fractal dimension (FD) analysis in Solea solea fed different experimental diets. 

Results: Histomorphological evaluation of all intestinal segments did not show signs of degeneration or 
inflammation. Cell proliferation index and FD were significantly reduced with a diet high in mussel meal 
(MM; p = 0.0034 and p = 0.01 063, respectively), while apoptotic index did not show any significant difference for 
the same comparison (p = 0.3859). Linear regression analysis between apoptotic index (independent variable) 
and FD (dependent variable) showed a statistically significant inverse relationship (p = 0.002528). Linear regression 
analysis between cell proliferation index (independent variable) and FD (dependent variable) did not show any 
significant correlation (p = 0.1 31 582). 

Conclusions: The results demonstrated that diets containing increasing levels of mussel meal in substitution of 
fishmeal did not incite a hyperplastic response of the intestinal mucosa. The mussel meal, which is derived from 
molluscs, could mimic the characteristics of the sole's natural prey, being readily digestible, even without increasing the 
absorptive surface of intestinal mucosa. Interestingly, from this study emerged that FD could be used as a numeric 
indicator complementary to in situ quantification methods to measure intestinal trophism, in conjunction with 
functional parameters. 
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PCNA, TUNEL, Fractal dimension 



Background 

The evaluation of intestinal trophism, mainly the mu- 
cosal layer, is an important issue in various conditions 
associated with injury, atrophy, recovery, and healing 
of the gut. The integrity of the mucosal barrier and the 
speed of either proliferation or apoptosis are other 
common variables of relevance for many investigative 
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studies [1]. To make these investigations objective and 
measurable, cellular turnover or morphometric param- 
eters such as villus height, crypt depth, and enterocyte 
height are frequently used [1]. Focusing on fish, histo- 
logical analysis of the digestive system is considered a 
good indicator of the nutritional status of health [2], 
The intestine and liver are the most important organs 
in digestion and the absorption of nutrients from food 
and, therefore, the monitoring of these organs is con- 
sidered necessary. Intestinal villi of mammals' and birds' 
digestive systems fit the Euclidean shape (a cylinder-like 
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structure) and allow the employment of the classical mor- 
phometric parameters. The morphological peculiarities of 
the fish intestine that distinguish it from the mammalian 
intestine, including the lack of distinct crypt compart- 
ments and the presence of longitudinal folding of the 
mucosa instead of villi [3], do not permit the use of this 
classical approach. 

The aim of the present study was to evaluate the kin- 
etics of the proliferation and apoptosis of enterocytes by 
immunohistochemistry and subsequently, to assess the 
complexity of intestinal mucosa by fractal dimension 
analysis in Solea solea fed different experimental diets. 

Methods 

The experiment was evaluated and approved by the 
Ethical-Scientific Committee for Animal Experimenta- 
tion of the University of Bologna, in accordance with 
the European Community Council directive (86/609/ 
ECC). Fish were obtained from the final sampling of an 
ongrowing experiment carried out at the Laboratory of 
Aquaculture, Department of Veterinary Medical Sciences, 
University of Bologna, Italy. Common sole {Solea solea) 
juveniles used for this experiment were obtained from 
natural spawning of the broodstock adapted to captiv- 
ity at the above-mentioned laboratory. From the larval 
stage, fish were cultured in a 3,000-L circular tank. 
Larvae were reared using the standard protocol adopted 
[4,5] and post-larvae and juveniles were fed a commercial 
diet (Biomar, Denmark, crude protein 63%, total lipids 
14%). Before the initiation of the trial, 840 sole juveniles 
were distributed into twelve 500-L square tanks (bottom 
surface: 0.64 m2) (70 individuals tank-1) in a completely 
random manner, to acclimatize them to experimental con- 
ditions for 41 days. During this period, fish were fed by 
automatic feeders for 12 h per day with a mixture of the 
four experimental diets. The tanks were provided with 
natural seawater and were connected to a unique closed 
recirculation system consisting of a mechanical sand filter 
(0.4 m3 of silica sand, 0.4-0.8 mm, Astral pool PTK 1200 
model filter, Servaqua S.A. Barsareny, Spain), an ultraviolet 
light (PE 25 mj/cm2: 16 m3 h-1, Blaufish, Barcelona, Spain) 
and a biofilter (PTK 1200, Astral Pool, Servaqua S.A. 
Barsareny, Spain). The water exchange rate per tank was 
100% every 2 h, whereas the overall water renewal of the 
system was 5% daily. Temperature was maintained constant 
at 20 ± 1°C throughout the experiment. The photoperiod 
was held constant at a 12-h day-length via artificial light 
(200 lx at the water surface, Delta Ohm luxmeter HD-9221, 
Delta-Ohm, Padua, Italy). The oxygen level was kept 
constant (7.5 ± 1.0 mgL-1) by a liquid oxygen system 
connected to software (B&G Sinergia snc, Chioggia, Italy). 
Ammonium (total ammonium nitrogen 0.5-0.1 mg L-l), 
nitrite (NO2-<0.2 mg L-l) and nitrate (N03- < 50 mg 
L-l) were determined spectrophotometrically, daily 



(Spectroquant Nova 60, Merk, Lab business, Darmstadt, 
Germany) at 12:00. At the same time, pH (7.8-8.2) and 
salinity (20 g L-l) were also determined. Each experimen- 
tal diet was fed to triplicate groups (13.1 ± 2.3 g initial 
mean body weight), assigned in a completely random 
manner, over 91 days. Fish were hand-fed twice-daily for 
1 h at a time (at 9:00 and 16:00) for 6 days a week and 
once on Sundays (9:00), to apparent satiation with four 
isoproteic (53%) and isolipidic (11%) palletized diets con- 
taining graded levels of mussel meal (0%, 25%, 50%, and 
75%, named respectively MM0, MM25, MM50, and 
MM75) to replace fish meal (FM). The meal distribution 
was made following a different tank order daily. To reduce 
the feeding hierarchies and to allow the access to the feed 
for all animals, the feed has been distributed evenly over 
the surface area. At the end of the trial, the average fish 
weight ± standard deviation (n = 3 tanks, 70 animals in 
each tank) was 32.4 ± 3.4, 42.4 ± 0.9, 46.0 ± 1.3 and 46.0 ± 
3.0 g in the groups fed MM0, MM25, MM50 and MM75, 
respectively. A linear regression analysis to analyze final 
weight against the MM dietary inclusion produced a R2 
value of 0.71 (p < 0.001). 

At the end of the trial, three fish per tank (total of 36 
animals) were sampled for gut histology. After euthan- 
asia with a lethal dose of 2-phenoxyethanol, the gut was 
removed and the intestine was divided into three seg- 
ments (anterior, intermediate, and posterior); from each 
segment, a 5 mm-long piece was sectioned and fixed in 
10% buffered formalin. Samples were then processed for 
routine histology to obtain a transversal section, which 
was stained with haematoxylin and eosin (H&E). Sec- 
tions were evaluated under a light microscope (Nikon 
Eclipse 80i, Nikon Corporation, Japan) for degenerative 
and inflammatory changes. 

Immunohistochemistry 

Anti-proliferating cell nuclear antigen (PCNA) antibody 
was employed for detection of the proliferation rate of 
enterocytes. Sections were deparaffinised for 30 min in 
Solvent Plus (Carlo Erba, Reagents S.r.L, Italy) and hy- 
drated in a graded series of alcohols. Subsequently, en- 
dogenous peroxidase activity was blocked with 3% 
hydrogen peroxide in distilled water for 30 min at room 
temperature. Sections were then treated in citrate buffer 
(pH 6.0) for 10 min in a microwave oven (750 W) for 
antigen retrieval. Nonspecific binding was blocked by 
incubating the sections with 5% normal goat serum and 
1% bovine serum albumin (BSA) in PBS buffer (blocking 
solution) for 1 h at room temperature and then sections 
were incubated overnight at 4°C with mouse monoclonal 
PCNA antibody (clone PC10, Cell Signaling Technology, 
Danvers, MA) diluted 1:4000 in blocking solution. A 
streptavidin-biotin-peroxidase polymeric complex (Super- 
Picture kit peroxidise, Zymed® Lab, San Francisco, USA) 
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was used to visualise antibody-antigen interaction and 
applied on the sections for 10 min at room temperature. 
The immunohistochemical reaction was developed using a 
3,3'-diaminobenzidine (DAB) solution (Sigma- Aldrich 
Corporate, St. Louis, MO). The slides were counter- 
stained with haematoxylin, dehydrated through alcohol, 
and cleared in xylene before mounting. 

Terminal deoxynucleotidyl transferase dUTP nick end 
labelling (TUNEL) assay 

The TUNEL method was used to detect the apoptosis 
rate of enterocytes. The ApopTag® Peroxidase In Situ 
Apoptosis Detection Kit (Millipore Corporation, Billerica, 
MA, USA) was employed. Sections were deparaffinised for 
30 min in Solvent Plus (Carlo Erba Reagents S.r.L, Italy) 
and hydrated in a graded series of alcohols. Treatment 
with proteinase K (20 ug/ml) in PBS for 15 min was per- 
formed. Subsequently, endogenous peroxidase activity 
was blocked with 3% hydrogen peroxide in distilled 
water for 5 min at room temperature. Equilibration 
Buffer was then applied on the sections for 3 min and 
finally, they were incubated with Working Strength TdT 
Enzyme for 1 h at 37°C in a humid chamber. After 
washing with Stop/Wash Buffer for 10 min, the second- 
ary antibody anti-digoxigenin peroxidase conjugate was 
applied on the sections for 30 min at room temperature. 
After washing, the sections were subjected to the same 
protocol as used for PCNA. 

Image analysis and cell counting 

From each stained section, photomicrographs of five 
fields per gastrointestinal tract at 20x magnification 
were captured using a Nikon Digital Sight SD-MS camera 
(Nikon Corporation, Japan) connected to an optical 
microscope. Images were then processed using ImageJ 
1.46 software, which is freely downloadable [6]. Acquired 
photographs of PCNA- and TUNEL-stained slides were in 
24 bit jpg format, with image resolution of 2560x1920 
pixels. 

PCNA-positive nuclei of enterocytes were counted 
using the plugin color segmentation [7]. The original 
image is divided into colour channels, which are manu- 
ally chosen by the operator as DAB-stained positive 
nuclei, haematoxylin-stained negative nuclei, cytoplasm 
of enterocytes, and white background. The software 
automatically calculates the percentage of area occupied 
by each chosen colour and, specifically, the percentage 
of area occupied by PCNA-positive (DAB-stained) 
nuclei. 

In order to eliminate false positivity of non-epithelial 
cells (erythrocytes, macrophages, lymphocytes, and 
fibroblasts), the lamina propria had been manually 
erased using the software Adobe Photoshop CS5 (Adobe 
Systems Inc, San Jose, CA). For apoptotic cell count, the 



following findings were considered to represent apoptosis 
[8,9]: (a) marked condensation of chromatin and cyto- 
plasm (apoptotic cells); (b) cytoplasmic fragments with 
or without condensed chromatin (apoptotic bodies); 
and (c) intra- and extracellular chromatin (apoptotic 
micronuclei). TUNEL-positive cells were counted manu- 
ally in order to eliminate the nonspecific signals due to the 
application of this method on intestinal sections [10-12], 
and then the plugin color segmentation was employed only 
to calculate the percentage of area occupied by the other 
colours (haematoxylin-stained negative nuclei, cytoplasm 
of enterocytes, and white background). 

All data obtained were used to calculate indexes as 
follows: 

n _ T1 . , % area PCNA-positive nuclei <#w% 

PCNA index = — — — * 100 

% area PCNA-positive nuclei + % area negative nuclei 

n° apoptotic nuclei/apoptotic bodies 

Apoptotic index = : — — 7 . —. ■ , ,. * 100 

n° negative nuclei per held + n° apoptotic nuclei/apoptotic bodies 

Where the n° of negative nuclei per field is equal to: 

total area of the field 20x * % negative nuclei 
mean area of a nucleus 

Image analysis and fractal dimension (FD) calculation 

Eighty photomicrographs of H&E-stained slides of intes- 
tinal tracts were obtained in blinded fashion at 4x mag- 
nification using a Nikon Digital Sight SD-MS camera 
connected to an optical microscope. Images were then 
processed using ImageJ 1.46 software: manual segmenta- 
tion and thresholding were performed blindly by a sec- 
ond operator to extract the one-pixel wide outline of the 
intestinal mucosal interface (outline) (Additional file 1). 
The FD calculation was performed with the Box Count- 
ing method using the specific command in ImageJ and 
setting boxes of 2, 3, 4, 6, 8, 1216, 32, 64, 128, 256, 512, 
1024, and 2048 pixels. Image resolution was 2560x1920 
pixels and format was 24 bit TIF for the original ac- 
quired images of H&E-stained slides; image resolution 
was 2560x1920 pixels and format was 8 bit Gif for the 
acquired photographs of outlines. 

Statistical analysis 

All data obtained were analysed by the software Statis- 
tica 8 (StatSoft Inc., Tulsa, Oklahoma). Normal distribu- 
tion was tested by Shapiro-Wilk. If normality criteria 
were met, one way ANOVA with post-hoc Fisher LSD 
test were run; if normality criteria were not met, 
Kruskal-Wallis multiple comparison were run. Signifi- 
cance was set as p < 0.05. 

Results 

Histomorphological evaluation of all intestinal segments 
did not show signs of degeneration or inflammation. 



Sirri et al. BMC Veterinary Research 2014, 10:148 
http://www.biomedcentral.eom/1 746-61 48/1 0/1 48 



Page 4 of 7 



PCNA labelling revealed a brown, intense, and homoge- 
neous nuclear staining of enterocytes mainly in the basal 
area and along the intestinal folds (Figure 1A). TUNEL 
labelling revealed a brown, intense staining of scattered 
single condensed nuclei or nuclear fragments (Figure IB). 

PCNA index was not normally distributed (n = 108; SW 
W = 0.96290; p = 0.00418). PCNA index in enterocytes 
showed a significant decrease in parallel with the increase 
in MM (n=108; H = 13.64672; p = 0.0034) (Figure 2). 
With respect to the intestinal tract, PCNA index did 
not show significant differences (n = 108; H = 4.298455; 
p = 0.1166). Apoptotic index was not normally distrib- 
uted (n = 108; SW W = 0.45922; p = 0.00000) and did 
not show significant differences with respect to the diets 
(n = 108; H = 3.037620; p = 0.3859). With respect to the 
intestinal tract, apoptotic index was significantly higher 
in the posterior segment compared with the anterior 
and intermediate segments (n = 108; H = 21.21437; 
p = 0.0000). FD was normally distributed (n = 80; Shapiro- 
Wilk W = 0.97847; p = 0.19599) and showed a significant 
decrease with a diet that was high in MM (n = 80; F 
(3.76) = 3.9997; p = 0.01063) (Figure 3). With respect to 
the intestinal tract, FD did not show significant 
differences (n = 80; current effect: F(2.77) = 1.3182; p = 
0.27359). Linear regression analysis between apoptotic 
index (independent variable) and FD (dependent variable) 
showed a statistically significant inverse relationship 
(Beta = -0.337225; p = 0.002528). Linear regression ana- 
lysis between PCNA index (independent variable) and FD 
(dependent variable) did not show a significant correlation 
(Beta = 0.164612; p = 0.131582). 

Discussion 

In the present study, the histomorphological evaluation 
of all three intestinal segments did not show signs of de- 
generation or inflammation, suggesting that MM could be 
considered a valid substitute for FM in the diet of sole. 
Immunohistochemistry with PCNA antibody and the 



TUNEL method have been widely used in fish species to 
demonstrate the renewal of enterocytes in studies ranging 
from experimental diet to environmental and toxicological 
assessments [13-21]. Focusing on the TUNEL method for 
apoptosis assessment, it has been widely demonstrated 
that it can produce an intrinsic nonspecific labelling. This 
does not represent a background staining but a binding 
with DNA fragments produced by endogenous cyto- 
plasmic endonuclease that is physiologically present, or 
by antigen retrieval pretreatment [10-12,22,23]. In the 
present study, the non-specific TUNEL labelling was cyto- 
plasmic, patchy, and hazy, with intense peripheral mem- 
branous positivity. Despite these apparent limitations, the 
TUNEL method coupled to H&E-stained sections enabled 
us to univocally identify cells undergoing apoptosis [12]. 

PCNA index showed a decrease in enterocytes prolif- 
eration with a diet with major inclusion of MM (MM75) 
with respect to the control diet (MM0), supporting the 
hypothesis of an hyperplastic adaptive response induced 
by a diet that deviates from the feeding habits of this spe- 
cies (a benthonic organism). On the other hand, apoptotic 
index was not influenced by the diet, probably because 
dietetic changes were not sufficient to induce a pro- 
apoptotic response compared with xenobiotics or hor- 
monal stimuli [13,14]. 

Moreover, the results obtained demonstrated a rela- 
tionship between reduced enterocytes proliferation and 
lessened complexity of mucosal folding, expressed by a 
lower FD, in the diet with the highest MM content. The 
biological interpretation of this result can be that diets 
containing high MM (MM75) are more similar to the 
natural feeding of sole. On the contrary, the control diet 
(MM0), which is usually administered in this reared 
species, is probably far from an optimal formulation; this 
could have triggered an hyperplastic adaptive response 
of the intestinal mucosa, supported by a higher prolifera- 
tion index, with the aim of increasing the absorptive sur- 
face. In linear regression analysis, an inverse and significant 



A 

£ 


B 


y 









Figure 1 Immunohistochemistry with anti-PCNA antibody and TUNEL assay. (A) PCNA-positive nuclei of enterocytes are located mainly in 
the basal area and along the intestinal folds. (B) TUNEL-positive apoptotic cells (arrows) and apoptotic bodies (arrow heads) are located at the 
apex of intestinal folds. The intrinsic nonspecific binding is also evident (asterisks) (Bars = 25 urn). 
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Figure 2 Diet effect on PCNA index. Pooled data of PCNA index (9 anterior + 9 intermediate + 9 posterior intestinal tracts = 27 for each group) 
for each diet displayed PCNA index of MM75 group significantly lower than MM50 and MMO (Post-hoc Kruskal-Wallis Multiple Comparison n = 108; 
H = 13.64672, significance p < 0.05) (letters mark significant differences in pairwise comparison). 



correlation between FD and apoptosis index emerged, 
confirming that a higher apoptosis rate corresponded 
to more simple mucosal folding even if there was no 
significant difference in apoptosis index in relation to 
the diet. 



1,38 



This study demonstrated the usefulness of FD calcula- 
tion as an alternative method to immunohistochemistry 
for the detection of intestinal trophism; FD permits a 
more rapid assessment of histological sections, reducing 
the costs of reagents and time. 
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Figure 3 Diet effect on intestinal fractal dimension. Pooled data of fractal dimension (respectively, MMO = 8 anterior + 7 intermediate + 8 posterior 
intestinal tracts = 23; MM25 = 4 anterior + 5 intermediate + 6 posterior intestinal tracts = 1 5; MM 50 = 6 anterior + 8 intermediate + 7 posterior 
intestinal tracts = 21 and MM75 = 6 anterior + 8 intermediate + 7 posterior intestinal = 21) for each diet displayed an increase of FD with the 
reduction of mussel meal (MM) diet content (Anova Post hoc Fisher LSD test n = 80; mean ± 1.96*SE, significance p < 0.05) (letters mark significant 
differences in pairwise comparison). 
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The FD calculation is able to transduce shape com- 
plexity in analytical quantitative data, harmonising the 
gap between structural features (typically qualitative) 
and functional quantitative measures. The fractal theory 
has been widely used in physics, biology, and medical 
fields and particularly in histology and histopathology. 
The FD of a pathological image can be interpreted as a 
snapshot of an evolving process; it is usually fractional 
and increases with the objects complexity [24,25]. In 
Euclidean shapes, a straight line has FD 1, whereas a 
fractal curve has a value of between 1 and 2 [26-28]. 
What the FD parameter provides is a statistical estimate 
of how an object fills a space across a broad range of 
scales; it is a useful parameter for the characterisation 
of complex, irregular structures, the analysis of which, 
when examined mathematically, denotes figures that 
resemble themselves when examined on different size 
scales (self-similarity) [29]. Focusing on histopathology, 
computing FD captures the architectural complexity 
of the histology specimen [27,30]. To the best of our 
knowledge, there have been no literature reports inves- 
tigating FD calculation of the intestine. The originality of 
the proposed method resides in the substitution or inte- 
gration of the conventional method for the detection of 
cell kinetics with a method of assessment of the complex- 
ity of the mucosal lining, thus intestinal trophism. 

Conclusions 

The results demonstrated that diets containing increas- 
ing levels of mussel meal in substitution of fishmeal did 
not incite a hyperplastic response of the intestinal mu- 
cosa. Fishes are not commonly present in the natural 
diet of sole, whose natural trophic profile is mainly 
composed of polychaetes and molluscs, as well as crus- 
taceans. The mussel meal, which is derived from mol- 
luscs, could mimic the characteristics of the soles 
natural prey, being readily digestible, even without in- 
creasing the absorptive surface of intestinal mucosa. 
Interestingly, from this study emerged that FD could be 
used as a numeric indicator complementary to in situ 
quantification methods to measure intestinal trophism, in 
conjunction with functional parameters. 

Additional file 



Additional file 1: Example of image processing and Box Counting 
FD calculation of an intestinal segment. 
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